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实数 编码 量子 共生 演算 法 及 其 在 云 任务 调度 中 的 应 用 
李 昆 仑 ， 关 立 伟 


(河北 大 学 电子 信息 工程 学 院 , 河北 保定 071000) 


摘 要 : 针对 共生 演算 法 收敛 慢 和 易 陷入 局 部 最 优 的 问题 ， 结 合 量 子 遗 传 算 法 理论 ， 提 出 一 种 实数 编码 的 量子 共生 演 
算法 (real-coded quantum symbiotic organisms search，RQSOS)。 首 先 依据 三 角 模糊 数 提出 差异 度 概念 ， 并 依 此 构造 一 个 
以 自 变 量 向 量 的 分 量 和 一 对 概率 幅 为 等 位 基因 的 三 倍 染 色 体 ， 使 一 条 染色 体 携带 更 多 信息 并 增强 解 的 多 样 性 ; 然后 提 
出 一 种 基于 阿 基 米 德 螺旋 线 的 探索 学 习 模 式 ， 加 强 对 解 空间 的 探索 精度 ; 最 后 使 用 共生 演算 法 更 新 差异 度 值 并 依据 差 
异 度 值 对 种 群 进 和 和 变异 操作 ， 促 使 整个 种 群 快速 向 最 优 方向 进化 且 减 小 了 陷入 局 部 最 优 的 概率 。 利 用 数值 优化 
问题 和 云 任务 调度 问题 对 算法 进行 验证 ， 仿 真 结果 表明 ，RQSOS 工法 在 收 伊 速度 和 寻 优 能 力 上 均 有 明显 提升 ， 是 一 种 
可 行 有 效 的 算法 。 
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Real-coded quantum SOS algorithm and its application in cloud task scheduling 


Li Kunlun, Guan Liwei 


(Electronic Information Engineering College Hebei University, Baoding 071000, Chinaa) 


Abstract: In order to solve the problem that symbiotic organisms search algorithm converge slowly and easy to fall into the 
local optimum, combining quantum genetic algorithm theory, this paper proposed a real-coded quantum symbiotic organisms 
search algorithm(RQSOS) . First, this paper presented the concept of the difference degree based on the principle of triangular 
fuzzy number, and constructed a variable component vector and a pair of probability amplitude of a allele in a chromosome that 
can catries more information and enhance the diversity of the solutions. Then the mode of rotary learning based on the 
Archimedes spiral was proposed, which strengthen the exploration ability of the solution space. Finally it is updated the 
difference degree based on SOS, and the population learning and mutation operations are carried out based on the value of the 
difference degree which can make the whole population evolution rapidly towards the optimal direction and reduce the 
probability of falling into local optimum. It is verified by numerical optimization and cloud task scheduling problem, and the 
simulation results show that the RQSOS algorithm can significantly improve the convergence speed and optimization ability, 
which is a feasible and effective algorithm. 
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scheduling 
0 引言 programming, QEP)、 量 子粒 子 群 算法 (quantum particle swarm 
一 optimization，QPSO) 和 量子 蚁 群 算 法 (quantum ant colony 


量子 遗传 算法 (genetic quantum algorithm,QGA)JID 是 algorithm, QACO) 等 B-5]。 

Narayanan 和 Moore 于 1996 年 提出 的 智能 优化 算法 , 首次 将 量 量子 智能 算法 能 为 排序 问题 、 背 包 问 题 、 路 径 优化 问题 和 
子 计算 相关 理论 与 遗传 算法 结合 ， 提 升 了 算法 性 能 。 随 后 Han 任务 调度 等 问题 找到 更 优 解决 方案 [510。 但 近 些 年 多 位 学 者 证 
等 人 中 将 其 扩展 为 量子 ; AN evolutionary 明 并 非 所 有 量子 智能 算法 都 适用 于 解决 连续 优化 问题 和 多 参数 
algorithm，QEA)， 促 进 了 量子 计算 理论 与 智能 算法 的 结合 ， 逐 优化 问题 ， 且 对 于 高 精度 高 维度 计算 问题 ， 通 过 对 量子 位 进行 
渐 形 成 一 个 新 的 研究 方向 --- 量 子 智能 算法 。 到 目前 为 止 ， 具 有 观测 而 得 到 的 二 进 制 编码 方式 很 容易 引起 染色 体 长 度 灾 难 ， 导 
代表 性 的 量子 智能 算法 有 量子 进化 规划 (quantum evolutionary ”和 致 求解 的 精度 降低 ,算法 进化 效率 下 降 0591。 为 克服 这 些 缺 点 ， 
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以 实数 编码 为 主要 思想 的 量子 智能 算法 受到 学 者 们 的 推 尝 。 晶 
于 实数 编码 量子 智能 算法 使 用 量子 比特 编码 染色 体 ， 利 用 量 了 
门 更 新 种 群 ， 其 通过 构造 量子 位 到 实数 部 的 映射 关系 ， 克 月 
传统 量子 智能 算法 解码 复杂 度 高 、 测 量 操 作 易 产生 退化 和 
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2 ”实数 编码 量子 共生 演算 法 


实数 编码 量子 共生 演算 法 是 一 种 融合 量子 遗传 算法 思想 和 
共生 演算 法 优点 的 智能 算法 。 该 算法 为 解决 共生 演算 法 前 期 进 


数 灾 ” 等 缺点 。 但 量子 位 到 实数 的 映射 关系 一 般 为 缩放 体制 ， 
解码 精度 不 能 得 到 保证 ， 故 这 些 算法 依然 存在 求解 精度 低 、 收 
敛 慢 和 易 陷入 局 部 最 优 的 问题 。 

为 了 克服 上 述 问 题 ， 结 合 量子 遗传 算法 和 共生 演算 法 ， 提 
出 一 种 实数 编码 的 量子 共生 演算 法 ， 利 用 数值 优化 问题 和 云 任 
务 调度 问题 对 算法 性 能 进行 验证 。 实 验 结果 表明 ， 该 方案 在 提 
升 共生 演算 法 性 能 的 同时 ， 成 功 避 免 了 上 述 缺 陷 。 


1 ”经 典 量子 遗传 算法 和 共生 演算 法 


1.1 经 典 量子 遗传 算法 
量子 遗传 算法 染色 体 用 量子 位 表示 , 一 个 量子 位 由 [wb8]7 组 
成 ,6 代表 状态 0 和 |1) 的 概率 幅 , 满足 归 一 全 和 件 |xjP2+IO2=1。 


化 速度 慢 和 后 期 种 群 多 样 性 丢失 的 问题 ， 首 先 改进 量子 比特 编 
码 形式 ， 构 造 一 个 具有 关联 特性 的 三 倍 染色 体 ， 然 后 通过 监视 
染色 体 信息 确定 种 群 进化 速度 和 多 样 性 情况 ， 并 依据 监视 结果 
适时 加 强 种 群 进 化 速度 并 补充 种 群 多 样 性 ， 最 后 利用 共生 演算 
法 更 新 染色 体 ， 实 现 问 题 的 优化 求解 。 
2.1 染色 体 编码 
本 文 提出 一 种 具 


关联 特性 的 编码 形式 ， 表 示 为 


一 个 量子 染色 体 可 表征 解 空 间 中 任意 解 的 亚 加 态 ， 有 具体 表示 如 


|. 和 Cn dl 
Pe ) 


量子 遗传 算法 依据 量子 门 进行 演化 ， 促 使 量子 种 群 趋向 最 
常用 的 量子 门 有 旋转 门 、 异 或 门 和 Hadamard 门 等 , 应 用 
较 广 的 是 如 下 的 量子 旋转 门 。 


a, cos@ -sinO | |a, 
| -ee | © 
其 中 : [wap 为 第 工 个 量子 位 ，[x28 下 为 旋转 之 后 量子 位 ，0 
为 旋转 角 ， 根 据 设 定 调整 策略 确定 其 大 小 。 旋 转 门 可 以 实现 状 
态 间 的 转换 ， 具 有 较 高 的 并 行 性 ， 保 证 算法 的 收敛 。 
1.2 共生 演算 法 
< 生 演算 法 
等 人 在 2014 年 提出 ， 
a) 共 生 过 程 。 
be ,= X,+rand(0,1)*(X,,, —Mutual_Vector* BF ) (3) 
Xj =X, ;+rand(01)* (x se — Mutual _Vector* BF,) (4) 


new 


(symbiotic organisms search, SOS ) {14 


包含 如 下 三 个 过 程 : 


Cheng 


Mutual _Vector =(X, + X,)/2 (5) 
其 中 :五代 表 第 i 个 位 置 ， 右 是 随机 与 名 相 作用 的 位 置 。BFj 和 
BF 代表 受益 水 平 取 值 为 1 或 2。 
b) 共 栖 过 程 。 
Xo = Xt rand(—1,1)*(X,, —X,) (0) 


其 中 :razd(-1,1) 为 [-1.1] 的 随机 数 ，Cer 芒 反映 一 种 受益 关系 ， 
节 提 供 一 些 受益 点 来 提升 名 的 存活 率 。 
co) 寄生 过 程 。 


Xinow = X， 
寄生 个 体 通 去 


过 随机 修改 不 同 维 数 中 的 数值 产生 ， 该 过 程 角 
加 强 种 群 多 样 性 ， 防 止 陷入 局 部 最 优 。 


x i 
cos(9) … cost, (8) 
sin(9) … sin(®,) 

其 中 : n 为 求解 问题 的 维 数 ，xi 为 自 变 量 向 量 x 的 分 量 ， 

%F=0.5*0*r，0 为 差异 度 值 ， 由 三 角 模 糊 数 产生 。 由 cos(J 控 制 

变异 信息 ，sin(9) 控 制 趋向 信息 。 

2.1.1 三 角 模糊 数 

依据 集 对 理论 5 , 三 角 模 糊 数 特征 可 由 二 元 联系 数 x+si 表 

示 ， 其 中 x 为 均值 ，s 为 方差， 计算 方法 如 下 : 

#3 + + (9) 


s=5e -二 + 人 -二 + 人 -二 (10) 
其 中 :x 为 三 角 模 糊 数 的 下 确 界 ,x" 为 三 角 模 糊 数 的 最 可 能 值 ， 
x" 为 三 角 模糊 数 的 上 确 界 。x 和 s 是 对 三 角 模糊 数 相对 确定 与 
不 确定 性 的 描述 ， 在 D-u 空间 可 表示 为 如 图 1 所 示 。 


不 确定 性 描述 


UCO=x X+Sl 
ee | 
0 卫生 性 反光 
X 
到 1 三 角 模糊 数 在 D-U 空间 的 解释 


其 中 9 表示 相互 作用 的 方向 ， 称 为 特征 参数 U(x) 的 幅 角 
方式 为 : 9 =arctans/x。 

2.1.2 种 群 差异 度 

2.1.1 节 分 析 可 知 ，0 值 越 小 (9E[0,m/2])， 确 定性 越 大 ， 
波动 性 越 小 ， 因 此 ， 种 群 差异 度 由 下 述 方式 得 出 : 


? 计算 


* -ak Xi 十 Xo 十 | (11) 


= [EA 一 x 十 人 一 加) 十 (wi 


其 中 : 代表 最 优 信息 的 均值 ， 本文 称 其 为 最 优 均值 ，x, 为 当 
前 维 数 最 小 值 ， 蕉 ,代表 当前 维 数 最 优 值 ， 世 ,代表 当前 维 数 


-x) (12) 
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最 大 值 ， 计 算 s 时 ，x” 为 当前 位 置 值 ， 
性 关系 。 故 差 异 度 值 定义 为 g=arctans'/x o 


以 下 举例 说 明 染 色 体 实数 部 与 差异 度 间 关 系 : 


某 一 维 多 样 性 
1 3 
2 3 
[3 Tw| 2 
4 3 
5 3 
前 期 后 期 
位 置 i 位 置 i 
图 2 染色 体 实 例 
通过 计算 得 到 : 


OF=[0.5236 0.4636 0.4405 0.4636 0.5236]’。 


均值 越 接 近 9 角度 值 越 小 ， 
种 群 后 期 多 样 性 较 小 ， 


如 图 


较 强 ， 设 3 为 当前 最 优 值 ， 如 图 2 左 所 示 。 


s=[1.7321 1.5000 1.4142 1.5000 1.7321]’; 


可 知 ， 当 前 值 与 最 优 


相反 则 越 大 ， 但 0 依然 较 大 。 
2 右 所 示 ， 通 过 


计算 得 


s=[0.4082 0.4082 0.5000 0.4082 0.4082];0=[0.1519 0.1519 


0.1853 0.1519 0.1519]’。 


可 知 ，9 依然 满足 与 最 优 


均值 越 接 近 值 


越 小 的 关系 , 算法 后 期 多 样 性 差 , 故 差 异 度 值 更 小 。 图 3 为 SOS 
算法 优化 Rosenbrock 函数 的 平均 差异 度 值 变 化 曲线 。 
allt 
os | | 
i 
Ml Wh IM 
省 [NN ap 
0 100 200 er 300 400 500 
图 3 SOS 算法 9 值 随 迭 代 次 数 变 化 曲线 
图 3 可 知 ， 种 群 后 期 多 样 性 已 经 严重 丢失 ， 算 法 性 能 下 
降 。 因 此 ， 通 过 对 0 值 的 监测 ， 实 时 了 解 种 群 多 样 性 ， 当 种 群 
多 样 性 过 大 时 ， 应 采取 措施 增加 个 体 优秀 基因 ， 可 加 快 算法 收 
敛 速度 。 在 种 群 多 样 性 丢失 时 及 时 采取 措施 可 减 小 算法 陷入 局 
部 最 优 的 概率 ， 加 强 算法 寻 优 能 力 。 


2.2 趋向 学 习 和 交叉 学 习 


2.2.1 基于 阿 基 米 德 螺旋 线 的 探索 学 习 模 式 


文献 [16] 和 [17] 证 


明了 反 向 学 习 和 旋转 学 习 的 有 效 性 。 


但 旋 


转角 度 的 设 定 影响 着 算法 探索 精度 ， 该 模式 依然 存在 探索 粗粮 
的 问题 。 为 进一步 加 强 算法 趋向 最 优 解 和 探索 能 力 ， 提 出 一 种 
基于 阿 基 米 德 螺 旋 线 的 探索 学 习 模 式 。 

如 图 4 所 示 ， 旋 转 模 式 分 为 两 种 ，a) 并 行 旋转 模式 ， 
所 加 半径 为 半径 绕 圆心 旋转 ， 半 径 与 圆 及 阿 基 米 德 线 相交 点 横 
坐标 即 为 并 行 学 习 点 ; nae 以 外 圆 半 径 为 半径 
圆心 旋转 ， 产 生 两 个 学 习 个 由 阿 基 米 德 线 内 部 最 小 圆 上 


> 
员 | | 


By 
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圆 上 移动 ， 这 样 一 个 在 外 


圆心 向 外 学 习 ， 可 加 1 


天 学 习 速率 。 具 体 过 程 如 


下 : 


半径 作 外 转 


以 C 为 中 心 (C=(a+b)/2，a,b 为 解 空间 上 下 界 )，(b-a)/2 为 


圆 并 以 式 (13 ) 在 设 定 空间 


Py 


半径 变化 范围 为 0~(b-a)/2。 


内 散发 阿 3 


基 米 德 螺 旋 线 ， 


f(n)=(a,b ,Ga #0 (13) 
y 4 
图 4 ”基于 阿 基 米 德 螺旋 线 的 探索 学 习 模 式 
tn) 代表 螺旋 线 的 产生 方式 ， 其 极 坐标 形式 为 -=ag+b 。 
学 习 过 程 为 对 当前 点 进行 一 定 角度 5 的 旋转 ， 旋 转 后 位 置 则 为 
学 习 后 的 点 ， 其 横 坐 标 为 
五 = [ 空 ]+yO cos(w+D) (14) 


其 中 :8 为 旋转 角 , a=arccos((2*z-q-bp)/(b-q))，flnn) 为 旋转 位 置 


I 


径 ; 外 围 圆 上 点 对 应 横 坐 标 为 
atb a—b 
a 7 ji 3 J oat) (15) 
基于 阿 基 米 德 线 的 两 种 学 习 模式 采取 不 同 的 方式 求 得 学 习 
点 ， 可 加 强 算 法 的 学 习 精 度 ， 提 高 搜索 效率 。 
2.2.2 趋向 学 习 
由 2.1 节 可 知 ,9 角 取 值 较 大 (Csin(9)>P), 种 群 多 样 性 较 好 ， 
但 与 当前 最 优 均值 的 综合 距离 相差 较 大 ， 代 表 最 优 信息 较 少 。 
为 加 强 算法 的 收敛 能 力 ， 加 快 算法 向 最 优 方 向 逼近 ， 本 节 引 入 
趋向 学 习 操 作 ， 使 个 体 在 一 个 扇形 空间 内 趋向 最 优 个 体 ， 最 大 
趋向 半径 为 R'"， 如 图 5 所 示 ，poev 为 当前 最 优 位 置 在 外 围 圆 半 
径 与 x 轴 夹 角 ，pi 为 当前 位 置 夹 角 。 
到 
趋向 后 位 置 了 
这 _ | | 
2 人 A Es th, 
| 
图 5 扇形 区 域内 的 趋向 操作 
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扇形 空间 范围 和 趋向 个 体 数量 可 根据 实际 情况 调整 ， 趋 向 


企 | 
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况 ， 具 体 过 程 如 下 : 


方式 可 由 式 (16)~(18) 表 示 。 a) 初 始 化 实数 部 种 群 和 最 大 旋转 圈 数 g， 旋 转 精度 a， 趋向 
A ED A 个 体 数 量 n:， 变 异 个 体 数 量 n2 和 变异 率 7 等 。 
AG=—sgn(A):AQ, -exp| -一 ”一 一 (16) 、 
VPme — Vf jin b) 通 过 模糊 三 角 函数 形式 ， 计 算 差 异 度 值 并 产生 三 倍 染 色 
和 编码 形式 。 
vp max = at) ci 体 的 编码 形式 
Ox Ox c) 用 SOS 算法 确定 个 体 在 解 空 间 中 的 位 置 x;， 更 新 差异 度 
er) ee x,) 0。 
Vfiwin = mi 本 18 
"| Ox/ 《2 dd) 学 习 操 作 : 
其 中 :V(x ) 为 评价 函数 00) 在 处 的 梯度 ，Ago 为 初始 旋转 (a) 趋向 学 习 。 由 式 (16)-(18) 产 生 趋 向 个 体 ， 若 新 个 体 适 
角 ， 具 体 大 小 视 问 题 确定 ，-signC4) 控 制 旋 转角 度 方向 ， 其 中 “应 度 优 于 当前 个 体 ， 则 以 新 个 体 代 蔡 当 前 个 体 。 
cos(D ee ss 、 i (b) 交叉 学 习 。 由 式 (19) 产 生 交 叉 学 习 个 体 ， 若 新 个 体 适 
nd 2 o 旋转 3 度 为 new™ A 二 和 A 9 其 中 4 为 
sin(A ) sil 6) de le 应 度 大 于 当前 个 体 ， 则 以 新 个 体 代替 当前 个 体 。 
加 速 系 数 ， 用 以 确定 角度 旋转 圈 数 ， 取 值 为 =d*2x，4 为 [0.9] 6) 变 异 操作 。 由 式 (20)~(21) 计 算 变异 信息 , 并 产生 变异 个 体 ， 
之 间 的 随机 整数 ，g 为 规定 旋转 最 大 圈 数 。 采 用 模式 2 进行 旋 ”随机 蔡 换 种 群 中 2 个体。 


转学 习 ， 取 ni 个 趋向 个 体 ， 若 趋向 个 体 为 有 效 个 体 则 接受 新 个 
体 ， 否 则 淘汰 。 
2.2.3 交 又 学 习 


为 进一步 加 强 算法 性 能 ,提出 一 种 交叉 学 习 模 式 ,具体 为 : 
随机 选取 M 个 竺 交叉 点 ， 采 用 式 (19) 决 定 是 否 向 当前 最 优 个 体 
进行 学 习 ， 学 习 有 效 则 接受 新 个 体 。 

j 员 沪 
We hesif (sin(9)) >P (19) 
Ns else 


其 中 : 9=0.5*0*x， 上 式 可 解释 为 : 0 值 较 大 ， 说 明 当 前 维度 信 
息 离 当 前 最 优 均 值 较 远 ， 应 大 概率 接受 当前 最 优 第 j 维 的 值 ， 
相反 则 小 概率 接受 当前 最 优 值 。 依 据 差 异 度 值 决定 是 否 学 习 ， 
有 助 于 种 群 学 习 的 成 功率 ， 加 快 种 群 进化 速度 。 
2.3 变异 

学 习 操作 引导 整个 种 群 向 当前 最 优 位 置 进化 ， 由 于 算法 的 
贪 禁 选 择 特性 ， 后 期 种 群 多 样 性 可 能 丢失 ， 算 法 易 陷 入 局 部 最 
优 解 。 基 于 此 ， 本 节 提 出 一 种 变异 操作 。 
当 0 值 较 小 ， 认 定 种 群 与 当前 最 优 具 有 极 强 的 相似 性 ， 多 


样 性 丢失 ， 在 此 时 进入 变异 操作 。 依 据 阿 基 米 德 探索 模式 ， 利 
用 如 下 公式 产生 变异 信息 : 
M, +sgn(e): f(r) ,lf (cos(9)) >77 oo0) 
Xi ， else 
m0 We e[0.5,1.5] 六 
else 


判断 是 否 达 到 终止 条 件 ， 达 到 则 进入 g)， 否 则 返 
多 算 法 终止 ， 记 录 最 优 数据 。 
由 上 述 步骤 可 知 , RQSOS 是 一 种 具有 反馈 更 新 机 制 的 算法 ， 
如 图 6 所 示 


回 c)。 


A 
o 


概率 学 习 ， 变 异 


Xi 


9 


SOS 


图 6 RQSOS 图 示 


染色 体 中 ，SOS 给 出 的 基因 取 值 具有 非 确定 性 ， 
概率 决定 具体 取 值 方式 和 大 小 。 故 RQSOS 模式 中 染色 体能 以 
一 种 县 加 态 来 描述 解 空间 ， 能 够 保持 解 的 多 样 性 。 


3 ”实验 与 分 析 


本 章 以 数值 优化 问题 和 云 工作 流 任务 调度 问题 为 例 ， 与 其 
他 算法 进行 对 比分 析 ， 验 证 本 文 算法 性 能 。 
3.1 数值 实验 及 分 析 

实验 采用 表 1 所 列 12 个 基本 测试 函数 验证 算法 性 能 ， 表 
中 马 为 取 值 范 围 ，B 为 最 优 值 ，D 为 函数 所 取 维 数 ，Threshold 
为 求解 精度 。 

实验 对 QPSO、 混 合 粒子 群 算法 PSO-GA)、 实 数 编码 量子 
差分 演化 算法 (RQDE) 和 SOS 的 优化 性 能 和 结果 进行 比 对 ， 本 
文 算法 参数 g=12,a=0.25,n1=5,ns=3,W=PFPe=rand, A00=0.047, 各 
算法 独立 运行 30 次 ,最 大 迭代 次 数 为 10000， 当 算法 达到 规定 


其 中 : sgn(*) 为 变异 方向 ，r=Q*pBw, pm=rand*n*1，a 为 精度 控制 

变量 , x 为 [1,g] 的 随机 整数 ， 为 变异 闵 值 。 

利用 模式 1 在 规定 空间 随机 勘测 出 nz 个 变异 个 体 并 随机 考 

换 种 群 中 m 个 个 体 〈 除 最 优 个 体外 )。 通 过 变异 操作 针对 性 的 

补充 种 群 多 样 性 ， 可 使 算法 不 易 陷 入 局 部 最 优 解 。 

2.4 算法 流程 
本 文 依据 SOS 算法 更 新 实数 部 进而 更 新 差异 度 0 值 , 然后 


精度 或 最 大 迭代 次 数 时 终止 ， 记 录 30 次 实验 平均 值 与 方差 
依据 实验 结果 记录 算法 运行 成 功 次 数 。 表 2 列 出 了 五 种 算法 的 
实验 结果 ， 其 中 w/W1 表示 本 文 算法 胜 过 、 持 平和 劣 于 对 比 算 法 
的 函数 个 数 ，M 为 均值 ， 严 代表 方差 。 
由 表 2 可 知 ， 五 种 算法 中 本 文 算 法 表现 最 优 ， 获 得 了 7 个 
则 试 函 数 的 全 局 最 优 解 ,11 个 测试 函数 最 优 解 , 而 基本 SOS 算 
法 只 获得 了 3 个 测试 函数 的 全 局 最 优 解 ,6 个 测试 函数 最 优 解 ， 


eh 


依据 差异 度 值 来 控制 该 维 信息 在 一 定 范围 内 的 学 习 与 变异 的 情 


且 在 12 个 测试 函数 上 RQSOS 性 能 均 优 于 等 于 基本 SOS 算 


1 A 会 、 人 


任务 调度 中 的 应 用 


@ | 四 | 
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法 ， 除 在 优化 f6 函数 时 本 文 得 到 优化 效果 微 差 于 QPSO 外 , 均 本 
优 于 其 他 算法 ， 证 明 本 文 规则 是 有 效 的 。 和 i a 
表 3 展示 了 五 种 算法 在 不 同 测试 函数 上 的 执行 成 功率 ， 可 i Oe 
以 看 出 RQDE、SOS 和 本 文 算法 RQSOS 在 执行 成 功率 上 优 于 fi Rosenbrock 30 [-30,30] 0 1.00e-001 
PSO-GA 和 QPSO 算法 , 且 RQSOS 算法 成 功率 最 高 。 说 明 SOS 下 Acldey 30 [32.32] oe 
算法 的 探索 能 力 强 、 鲁 棒 性 好 的 优点 ， 也 进一步 说 明了 本 文 提 和 ee 
出 算法 保持 了 SOS 算法 优点 , 克服 了 算法 容易 陷入 早熟 的 缺点 ， PR 
能 达到 更 好 的 优化 效果 。 fi Dixon-price 30 [-10,10] 0 1.00e-006 
图 7 和 8 展示 了 五 种 算法 在 优化 函数 户 和 上 万 的 收敛 曲线 对 f Zakharov 10 [5,10] 0 1.00e-006 
比 图 。 由 图 可 知 ， 在 不 同 函数 优化 的 问题 上 ， 不 同 算法 的 展现 a 
了 不 同 的 性 能 ， 本 文 算法 均 有 较 高 的 收敛 能 力 和 较 强 的 探索 能 人 i i 
力 .RQSOS 能 产生 优良 解 且 具 有 较 好 收敛 能 力 的 原因 为 采用 了 0 a 
学 习 模式 和 变异 模式 ， 学 习 模 式 加 快 种 群 向 最 优 解 方向 发 展 ， 和 ee 
而 变异 模式 保持 了 种 群 的 多 样 性 ， 故 本 文 算法 能 够 均衡 算法 的 a 
勘测 能 力 和 收敛 能 力 ， 加 强 算法 性 能 。 
亏 表 2 本 文 算法 与 其 他 四 种 算法 的 仿真 结果 比较 
< Funs. PSO-GA MEV QPSO MEV RQDE MEV SOS MEV RQSOS MiV 
fi 2.83e-00+4.01e-00 6.83e-25+1.98e-48 1.28e-18+2.16e-36 3.01e-134+3.63e-267 9.20e-169+1.02e-280 
ft 3.23et+02+1.13e+04 2.49et+01+t1.54e-01 7.69e-00+1.39e+01 2.60e+01+1.41e-00 1.71e-3+6.52e-06 
fs 1.99e+01+0.00e-00 1.79e-14+1.21e-29 3.64e-10+1.69e-20 8.89e-16+3.89e-62 8.89e-16+3.9e-62 
fa $5.76e-01+8.78e-02 3.12e-02+5.28e-04 1.09e-12+2.91e-26 0.00+0.00e-00 0.00+0.00e-00 
fs 9.95e-02+8.91e-02 2.21et+01+4.57et+01 9.95e-02+8.91e-02 9.95e-02+8.91e-02 0.00+0.00e-00 
fo 3.96e-00+1.61e-00 0+0.00e-00 8.10e-03+5.90e-04 3.96e-00+1.61e-00 1.09e-05+4.70e-06 
fy 1.20e+01+t1.31et+03 3.33e-61+9.45Se-121 3.98e-66+1.19e-131 4.27e-290+0.00e-00 0.00+0.00e-00 
fs -9.6171+1.20e-03 -6.8226+7.88et+01 -9.6396+3.87e-04 -9.6248+9.62e-04 -9.6600+1.4Se-07 
fo 8.1062+3.16e-30 5.5428+2.83e-00 8.1062+3.1Se-30 8.1062+3.16e-30 8.1062+0.00e-00 
fio 3.01e-02+0.00e-00 0+0.00e-00 0.00+0.00e-00 0.00+0.00e-00 0.00+0.00e-00 
fi 0.00+0.00e-00 0+0.00e-00 0.00+0.00e-00 0.00+0.00e-00 0.00+0.00e-00 
fi -186.7308+1.04e-08 -186.7283+S.84e-005 -186.7309+0.00e-00 -186.7309+0.00e-00 -186.7309+0.00e-00 
W/t/l 11/1/0 9/2/1 8/3/0 7/5/0 - 
表 3 五 种 算法 的 执行 成 功率 对 比 
二 PSO GA QPSO RQDE SOS RQSOS ee 
SR% 
fi 17 100 100 100 -100 i 
f 0 0 30 10 87 i , 、 
fs 0 100 100 100 100 - ' 
各 37 23 100 100 100 ob ee ee 
f 23 17 87 100 100 
下 100 而 17 i 图 7 Rosenbrock 收敛 曲线 对 比 图 
J 图 9 和 10 为 SOS 和 RQSOS 算法 在 优化 函数 户 和 户 时 种 
号 0 群 差 异 度 的 对 比 , 可 知 SOS 算法 优化 户 函 数 时 种 群 多 样 性 很 快 
丢失 ， 探 寻 最 优 解 的 能 力 下 降 ， 故 其 优化 精度 和 执行 成 功率 均 
3 很 低 ; 在 优化 访 时 ， 虽 多 样 性 可 保持 在 一 定 范围 ， 但 其 向 最 优 
2 靠近 的 速度 很 慢 ， 因 此 ， 在 规定 选 代 次 数 内 很 可 能 不 能 搜寻 到 
最 优 解 。 很 明显 ， 本 文 策略 在 保持 高 收敛 性 的 同时 能 够 将 种 群 
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多 样 性 控制 在 一 定 范围 ， 进 一 步 证 明了 算法 的 有 效 性 。 
. QPSO 
PSO-G, 
-3 | RQDE 
sos 
者 汉 RQSOS 
®@ -5 
BE 
二 了 
a 
| 
-9 和 
9 50 100 200 5 300 350 向 450 和 
lter 
图 8 Michalewicz10 收敛 曲线 对 比 图 


2 


Difference degree 
oo 0o 


Difference degree 
o oo oo 


区 | 


500 


1000 


lter 


1500 


-一 80s 
Rasos 


2000 


10 处 理 妇 函数 SOS 与 RQSOS 算法 6 值 对 比 


3.2 ”本 文 算法 在 云 工作 流 任 务 调度 中 的 应 用 


目前 云 环境 中 工作 流 调 


可 


n 为 任务 数 ， S 


nan 辐 的 SyES 表示 人 牺 


本 文 将 处 


里 成 本 描 i 


述 为 人 有 


图 11 


为 一 个 简单 DAG 图 ， 


如 果 任 


应 ninj 间 的 通信 和 量 。 
处 理 器 间 的 通信 人 花费， 处 于 相同 处 至 


是 0。 


里 ， 从 而 缩短 进程 的 总 完成 时 间 ， 
DAG(DAG={N,5}) 图 描述 ， 
5 个 边 的 外 
E 务 三 必须 在 各 
E 务 的 执行 时 间 和 
非 负 权重 
务 依赖 于 不 同 


度 算法 从 不 同 角度 对 工作 流 进行 处 
提升 系统 性 


能 。 工 作 流 任务 


其 中 N 


n 个 节点 集合 构成 ， 


合 构成 ， 代 表 任 务 间 


E 务 ni 


的 约束 条 件 。 
完成 之 后 才能 执行 。 

使 用 处 理 器 的 费用 。 
Cij 关联 Sj ES 边 对 
的 处 理 器 ， 则 需 计算 


器 间 任 务 的 通信 花费 认为 


有 的 DAG 


李 昆 仑 ， 等 


处 于 处 理 器 P; 上 的 任务 


Tan (ni, 


则 任 


其 中 : Wlm) 为 任 
务 ni 为 任务 nj 的 前 


GhinaX 


实数 编码 量子 共生 演算 me 


ni 的 最 早 的 开始 时 


间 Tsart 


太刀 为 入 口 任务 ; 
LT, sl 站 else; 


pj)= | 要 人 


务 ni 的 完成 时 间 为 
FT， (n ,)= 


务 Ni 计算 量 了 里， Wpj) 代 表 处 


Taalnis pj) + Wm) Wlp,) 


里 器 的 处 理 能 


合作 期 刊 


调度 中 的 应 用 


nap)) 定 


(22) 


(23) 
才 。 任 


前 继任 务 ninj 所 在 


务 ， Popy 分 别 为 任 


元 ， 若 订 则 COzazh=0， 故 户 到 达 记 的 时 间 为 


AT, (n, 


EST, ln 


其 中 : 


ng 为 处 理 器 上 最 后 完成 的 任 


n))=FT, (n,)+ Cln,n;)/ Clp,,p,) 
))= max fmax {AT, (n,n)) FT, (n,)) 


务 ， 在 相同 处 理 


器 上 


处 理 音 


(24) 


Q5) 


任务 执 


行 的 优先 级 关系 


式 (26) 确 定 。 


Pri(n,)=W( 


n;)/M, +, max, (Cln, 


n,)/M, +Pri( 


:BL(n,) 


不 由: Mp,Mi 为 处 


表 任 务 ni 的 后 继 点 集合 ， 

故 本 文 的 适应 
[0,1], x+p=1,s 为 调整 系数 , 7yj 为 占 
器 单位 时 间 的 使 用 价格 。 
5 权重 。 


整 时 间 和 花费 所 


理 器 处 理 和 


nm))+ 


(26) 


E 力 中 值 和 传输 能 力 中 值 succCo) 代 
BL(n7) 为 任务 ni 的 Bottom-Level。 


度 函 数 可 以 由 公式 (27)-(29) 表 示 , 其 中 
第 j 个 处 理 器 的 总 


局 op E 
时 间 ， 


可 根据 实际 应 


Fitness=Q:Time +B:e:SP 


为 验证 本 文 算法 在 处 理 云 工 作 流 调度 问题 时 
处 理 器 和 任务 数 变 化 
QGA 和 改进 量子 粒子 群 算法 (Improved Quantum Particle 
Optimization, IJQPSO) 进 行 对 比 。 实 验 中 ， 算 法 最 大 和 迭 
400 或 600 次 ，QGA 和 IQPSO 种 群 规模 为 60， 本 文 算法 种 群 
个 体 数 
节 相 同 。 


规模 为 30， 趋 向 
7， 其 他 参数 与 上 


Time = max (gs Ls (n, ) 


SP= > Pp 


的 有 效 
时 的 调度 问题 进行 模拟 仿真 ， 


目 为 3， 变异 个 体 数 


5， 圈 


本 次 实验 采 


实验 室 设计 的 云 任务 调度 仿真 系统 ， 


户 革 


制定 参数 ， 随 机 
系统 。 节 点 计算 量 控 


生成 10-100 个 节点 的 任务 和 4-10 个 
在 [50,200] 内 , 节点 间 通 信 


宇 吊 


力 
力 


取 值 范围 [30,10 
取 值 范围 [10,30 
如 图 12 和 


可 以 看 出 IQPSO 算法 寻 


内 ， 节 点 最 大 入 度 和 出 


8; 资源 
[110]， 


计量 控 人 
度 为 10， 图 最 大 深度 为 
0], 资 源 单位 时 间 使 用 价格 范 国 
]。 
13 为 当 资 源 数 固 定 为 8 时 的 仿 
得 的 完工 时 间 效 果 比 QGA 差 ， 


处 理 


器 价格 要 优 了 


F QGA， 而 本 文 算法 在 时 间 和 花费 上 


制 在 [5,30] 


结果 ， 


用 而 调 


Q7) 


(28) 


(29) 


性 ， 对 
分 别 与 
Swarm 


数 9 为 


依据 用 
资源 的 


计算 外 
传输 色 


由 
但 使 用 
均 优 于 
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IQPSO 和 QGA 算法 。 调 
能 产生 求 得 的 完工 时 间 与 
更 好 的 均衡 性 。 


度 算法 在 约束 条 件 不 同 的 情况 下 ， 可 
花费 不 均衡 的 现象 ， 而 本 文 算法 具有 


李 昆 仑 ， 等 


实数 编码 量子 共生 演算 法 


图 16 和 17 展示 了 任务 数 ， 


值 变化 的 情况 ， 


学 习 模式 和 变异 模式 是 有 效 
在 任务 数 增 加 时 ， 任 务 


结果 可 知 ， 
期 进化 效率 较 高 ， 但 IQPSO 算法 性 能 
保持 较 高 进化 效率 且 优 于 IQPSO 算法 , 说 明 本 文 算法 所 提出 的 


的 。 
处 理 


IQPSO 算法 和 RQSOS 算法 在 前 
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很 快 下 降 ，RQSOS 算法 


变 得 复杂 ， 算 法 想 寻 得 一 个 稳 


定点 均 需 要 一 定 和 迭代 次 数 , 人 
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图 17 可 看 出 , 本文 算法 能 够 很 
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快 的 找到 较 优 解 ， 在 外 
局 部 最 优 的 能 力 ， 但 是 效率 极 


章法 后 期 ，IQPSO 和 QGA 算法 均 有 跳出 


氏 ， 本 文 算 法 在 后 期 以 更 大 概率 
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任务 数 固定 时 ， 完 工时 间 和 最 小 花费 


结果 ， 可 看 出 在 资源 数 少 于 6 时 本 文 算 
法 得 到 的 完工 时 间 多 于 QGA 和 IQPSO 算法 ， 当 资源 数 超过 

时 间 较 少 ; 
均 优 于 QGA 和 QIPSO 算法 ， 故 在 任务 数 固 
本 文 算法 依然 具有 较 优 性 能 


对 于 最 低 花 费 ， 本 文 算法 所 得 结果 
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步 说 明了 本 文 算法 的 优越 性 
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探索 学 习 模 式 ， 
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通过 监测 差异 度 值 实时 了 解 种 群 的 多 样 
较 强 但 综合 距离 当 


距离 , 当 种 群 多 样 性 


采用 趋向 学 习 和 交叉 学 习 模式 加 
种 群 多 样 性 丢失 ， 为 防止 算法 陷 


加 强 了 算法 的 探测 能 
果 表 明 ， 本 文 算法 


论 与 共生 演算 法 ， 提 出 一 种 实数 编码 


E 演 算法。 首先 提出 差异 度 概 念 和 基于 阿 基 米 德 线 的 
可 使 算法 在 解 空间 快速 、 细 致 的 搜索 最 优 解 。 
性 和 与 当前 最 优 均 值 的 


前 最 优 均 值 距离 较 远 时 ， 
快 种 群 向 最 优 方向 的 进化 ; 当 
入 局 部 最 优 ， 引 入 变异 操作 ， 


力 。 数 值 实验 和 云 工作 流 任务 调度 实验 结 
了 效 的 提升 了 SOS 算法 性 能 ， 能够 得 到 更 优 


解 。 但 监测 信息 


是 以 反馈 方式 作 | 


于 算法 ， 该 方式 增加 了 原 算 


法 复杂 度 , 因 此 ,在 以 后 的 研究 工作 中 有 必要 进一步 分 析 模 型 ， 
加 强 本 文 提出 思想 与 原 算法 的 融合 程度 ， 减 小 算法 的 复杂 度 。 
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